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Abstract: 

Approximate mean-field equations of motion for the classical chiral field are developed 
within the linear a model by means of a Hartree factorization. Both the approximate and 
the unapproximated equations of motion are augmented with a Rayleigh cooling term to 
emulate a uniform expansion, thereby allowing the extraction of observables relevant to the 
detection of disoriented chiral condensates, specifically the pion power spectrum, the pion 
correlation function, and the distribution of the neutral pion fraction. While the mean- 
field dynamics appears to be sufficiently accurate to be practically useful, the results also 
illustrate some difficulties associated with identifying the phenomenon experimentally. 
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1 Introduction 



High-energy collisions of hadrons or nuclei may produce extended regions of space 
within which chiral symmetry is temporarily nearly restored and the subsequent non- 
equilibrium relaxation towards the normal vacuum may then produce disoriented 
chiral condensates, a coherent oscillation of the pion field that is manifested through 
anomalous pion multiplicity distributions |l|, 0, ^ |, |^ . In order to assess the prospects 
for such DCC phenomena to actually occur and be detectable above the background, 
it is necessary to perform extensive dynamical calculations. Most dynamical studies 
have been carried out within the framework of the linear a model 0, with the chiral 
degrees of freedom being represented by a classical field H, |T0|, [11], |T2|, |13|, |1^, 



151 , |T6|, 0. For recent reviews of the topic, see refs. [|^, |20|| . 



The purpose of this paper is two-fold. The first is to derive approximate equations 
of motions by means of the Hartree factorization technique, by proceeding in a manner 
analogous to what was done recently in the context of statistical equilibrium . 
The resulting mean-field description is practically simpler to treat and it provides 
a useful conceptual framework for understanding the dynamics. Additionally, we 
wish to employ the linear a model for the calculation of key observables relevant to 
the detection of disoriented chiral condensates, namely the power spectrum of the 
emerging free pions, their correlation function, and the associated distribution of the 
neutral pion fraction. For this purpose, we emulate a Bjorken-type expansion in D 
dimensions by augmenting the equations of motion with a simple Rayleigh cooling 
term. Starting from thermal equilibrium, we then follow the dynamical evolution of 
the chiral field until large times when the asymptotic regime of decoupling has been 
reached and the observables can be extracted. Comparisons between the results of the 
mean-field dynamics with those of the non-approximated equation of motion provide 
a basis for assessing the approximation which appears to be sufficiently accurate to 
be of practical utility. Moreover, the calculated observables are interesting in their 
own right, as they illustrate the key features of the DCC phenomenon and bring out 
some of the difficulties inherent in the experimental search for signals. 

The presentation is organized as follows. In sect. ^ we develop the mean-field 
equations of motion and, in addition, we indicate how the approximation can be 
utilized for statistical equilibrium. Subsequently, in sect. |^, we describe how field 
configurations corresponding to various physical scenarios of interest can be prepared 
and treated. We then turn to the discussion of the extracted observables, sect. ^ and 
finally make our concluding discussion in sect. ^ 
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2 The mean-field approximation 



The linear a model offers a convenient tool for exploring the non-equilibrium 
dynamics of DCC formation. In the usual approximation, it describes the chiral 
degrees of freedom by means of a real classical field which forms an 0(4) vector, 
4>{r) = (cr(r), 7r(r)), and is subject to a non-linear interaction, V = ^{(jP' — v'^Y — Ha. 
(We use units where % and c are unity, and 0^ = (/)o0 = cr2 + 7r-7r where o denotes 
the 0(4) scalar product.) The equation of motion for the chiral field 0(r,t) is then 

[U + X{<f - v^)]ci, = Ha, (1) 

where the three parameters. A, f , H are adjusted to ensure reproduction of the ground- 
state field strength 0vac = A (=92 MeV) and the free meson masses (=138 MeV) 
and m„ (=600 MeV). 

For the present study, we enclose the system in a cubic box with volume VL = Li^ 
and impose periodic boundary conditions. It is then possible to uniquely decompose 
the field, 

0(r,t) = ±(t) + 50(r,t) = E '^kW e^''-'' , (2) 

k 

V'(r,t) = ±{t) + 5xl^{r,t) = Y^i^^it) e'"^-"- . (3) 

k 

The underscored part is the spatial average over the volume of the box, (p =< (p > 
and "0 =< V' >• It "will be referred to as the "order parameter", while the remaining 
fluctuating part represents quasiparticle excitations relative to that constant back- 
ground field. The coefficients in the Fourier expansions shown on the right are simply 
given by spatial averages, (f)^^ =< <f)exp{—ik-r) > and t/^^^ =< -0 exp(— ik- r) >. The 
terms corresponding to = then represent the order parameter </)q = (p and its 
time derivative ipQ = ip. Moreover, since the fields are real, the Fourier coefficients 
satisfy symmetry relations, (j\ = (p_i^ and ipl^ = "jA-k- 

2.1 Equations of motion 

By inserting the decomposition <p = <p + Sep into the equation of motion (^ and 
subsequently taking the spatial average, it is possible to derive an equation for the 
order parameter, 

+ Mo ^ = Ha , (4) 
where the 0(4) tensor is given by 

M = A(< 02 > -v^)I + 2A < 6(f>6cf) > , (5) 

with I being the 0(4) unit tensor. We note that < 0^ >= 0g + < 50^ > where 0o is 
the magnitude of the order parameter, 0q = o 0. The equation (|D was first derived 
by Baym and Grinstein |^. Although exact, it is of little practical value unless the 
averages entering in the mass tensor (^) are provided. 
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Subtracting now the equation for the average field strength from the full equa- 
tion of motion (|l|), we obtain an approximate equation of motion for the field fiuctu- 
ations. It can be written on a simple Klein-Gordon form, 

□50 + M o 5(f) ^ , (6) 

where the quasi-particle mass tensor M is 

M = A(< 02 > -v^)I + 2X<(f>(f)> = M + . (7) 

The last relation follows because < (p(p > = (fxp + < 6(p6(p >■ In deriving this 
result, we have employed a Hartree linearization. Using j =0,1, 2, 3 to denote the 0(4) 
components of the chiral field, the approximation can be characterized as follows: 
Terms containing an even number of fiuctuations are replaced by their spatial average, 
S(j)j6(f)j' — 6(j)j6(f)j' > (which leads to a cancellation of the even terms) and the odd 
fiuctuation terms are contracted, 

6(j)j6(j)jr6(j)jrr < 6(f)j6(j)j/ > 6(f)j// + < S S > S(j)j + < S(f)'jS(j)j > S(f)j/ . (8) 

Generally, one could imagine an effective theory containing higher powers of the field 
(although it could not have arisen from a renormalizable quantum-field theory). The 
linearization method can be applied in this more complicated well by simply 

employing the gaussian approximation to reduce higher-order terms. 

Since the present treatment defines the order parameter as the average of the 
field amplitude over the entire volume of the box, it has no spatial dependence and 
the d'Alambert operator □ then reduces to df. However, it is possible to generalize 
the definition of the order parameter so that it does vary in space (see sect. |3.2|). The 



above mean-field equations can then still be employed and with a view towards this 
more realistic situation we leave the full d'Alambert operator in place. 

The above developments are in accordance with earlier work by Boyanovsky et 
al. ||13|. However, in so far as the dimensionality of the chiral space can be regarded 



as large, ^ 1, its inverse provides a useful expansion parameter, and it may 
then be natural to keep only the leading terms in 1/A^ [0, ||, [l^, |13|, |16|, |4 



The term 2A < S(f)S(f) > is then absent from the mass tensors iV£ and M given in 
eqs. (^ and (|^). This term is of relative order 2/A^, since the leading term in the 
effective masses contains A^ individual fluctuation terms < 5</)| >, and there are no 
higher-order terms. Since the actual dimensionality is not all that large, N=4, the 
inclusion of the 1/A^ terms produces a 50% increase in the contribution from the 
quasi-particle degrees of freedom to the effective masses. Changes in the fluctuations 
of such relatively large magnitude may have a signiflcant influence on the calculated 
results and we therefore prefer to keep all the terms. (It should be noted that the 
inclusion of the term 2A < 5050 > has no impact on the computational burden, 
since it merely modifles the magnitude of the appropriate coefficient.) As we shall 
see, the mean-fleld treatment with the 1/A^ term retained provides a practically useful 
approximation to the full equation of motion (0). 
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2.1.1 k representation 

The mean-field description derived above provides a very convenient means for ob- 
taining approximate solutions for tlie dynamical evolution of the chiral field. When 
reexpressed in the k representation by means of a Fourier transformation, the coupled 
equations of motion are as follows, 

k = : df(f)Q + Mo(f,Q = Ha , (9) 

k>0 : {dl + k^)(j)^ + Mo0j^ = 0, (10) 

where the mass tensors are given by 

M = A(^0^o0^ - v')I + 2A^0k0k = M + • (H) 

k k 

The above equations (^0) form a simple self-contained system that can be solved at 
only modest numerical effort. In sect. ^ we shall illustrate the utility of this convenient 
approximate description. 



2.2 Energy 

The total energy associated with a given dynamical state of the chiral field can be 
written as QE[ip, (f)], where </>] is the average energy density. This key quantity 
can be decomposed into terms having instructive physical interpretations, thereby 
enabling us to develop a simple approximate treatment of statistical equihbrium 

E[il^, 0] = < + 1 (v^)2 ^ A^^2 _ ^2)2 _Ha>= Eo + E^^ + SV . (12) 

Here Eq is the energy density arising if all the field fluctuations were put to zero, 

Eo ^ E[V = V^(r),0 = 0(r)] = + ^(00-^')' - ^0ocosxo, (13) 

where t/'q = i/j o i/j and the disorientation angle Xo measures the angle between the 
order parameter (p and the vacuum direction a, cosxo = 4> o cr. The first term in Eq 
is the bare kinetic energy density Ko and the remainder is the bare interaction energy 
density Vq. 

The second term in (|l^) is given by 

£^qp[i/',0] = -<6ipo6ip + V6(poV6(t) + 6(j) o M o 6(f) > , (14) 

where M is the quasiparticle mass tensor (0) derived in connection with the mean- 
field equations of motion above. Therefore fi-Eqp can be interpreted as the energy 
associated with the gas of independent quasiparticles described by the Klein-Gordon 
equation of motion d). 
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The last term in ([T2|) corrects for the fact that the interaction is non-hnear. It 
can be written as 

W[il^, 0] = ^<5<i)^> - ^ < 50^ >G + X±o<54) 5<i? > , (15) 

where < ■ >g denotes the evaluation of the average by means of the gaussian approx- 
imation, 

< >G = 5(t)j5(t)j X 6(j)j>6(f)j> > + 2 < 6(f)j5(j)j> >< 6(f)j>5(f)j >] 

jj' 

= < (50^ >2 + 2 tr(< > o < S<t)6(t) >) . (16) 

It is important to recognize that the above relations are exact, since they are merely 
reformulations of the expression (|I2|) for the energy of a given state. 

2.3 Statistical equilibrium 

We now briefly indicate how the mean-fleld approximation can be applied to statistical 



equilibrium, leading to the simple approximate treatment developed in ref. pi 
The natural starting point is the partition function. 



|l?[V',0]e-^^[^''^], (17) 



where the functional integral is over all physically possible conflgurations of the chiral 
field ('i/j(r'), (/)(r)). For the purpose of treating the statistical equilibrium, we now 
make a number of approximations. The first and most important approximation is to 
replace the state-dependent correction term 6V ( |T5|) by its thermal average, for each 
given value of the order parameter {\l),(f)). This is a good approximation when the 
system considered extends over a volume that is large on the scale of the correlation 
length, since the fiuctuations in one part of the system are then independent of those 
in other parts and so the spatial average < ■ > approaches the corresponding thermal 
average -< ■ 



The last term in (|1^) then vanishes by symmetry, since it contains an odd number 
of field fiuctuations. Moreover, we employ the gaussian approximation (^) for the 
evaluation of the first term in (p^). We can simplify the analysis by taking advantage 
of the fact that, for any given value of the order parameter (p, the tensor -< 5cf)5cf) >- 
is aligned with the order parameter and, furthermore, is invariant under rotations 
around that direction. It then follows that 

tr(^ 6^6(f) = ^ (50V = ^ 50j ^ + 3 ^ (50i ^ , (18) 

tr(^ 6cj)6(f) S(f)S(f) ^) = ^ 50j ^2 _^ 3 ^ ^2 ^ (^g) 

where (50|| denotes the component of 64> in the direction of the order parameter, 
(5011 = 6(p o <P and -< 505_ >- denotes the variance of the thermal field fiuctuations in 
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(any) one of the perpendicular directions. These approximations amount to a simple 
replacement, 

5V['iJ^,4>] ^ 5V{±,^) = -^^50Vg= (20) 

[3 ^ 6(1)1 + 6 ^ 50j 50i ^ + 15 ^ y^] . (21) 

Once 6V in the energy expression (|12D has been approximated in this manner, the 
partition function can be brought onto a convenient factorized form. 



e-T^^p , (22) 



where the functional integral over the states of the chiral field has been expressed in 
terms of the Fourier components introduced in eqs. (@-^. The prime indicates the 
omission of the zero mode having k = 0. Since the state-dependent quasiparticle 
energy E^p is a sum of contributions from each wave number k, 

^qp = ^E'[V'k°V'k + ^'0k°0k + 0k°^°0k] , (23) 

^ k 

the multiple integral over the fluctuations factorizes correspondingly. 

The expression for the partition function may then be reduced to an integral over 
the order parameter (and its time derivative), 

Zt = j <f±d'^±WT{±,^ , (24) 
with the associated statistical weight being given by 

1^t(V',</>) ~ Q-l^iEo+SV+'iE^^^-TST) ^ ^-q:{Ko+VT-TST) ^ g-§(ifo+FT)^ (25) 

The last two terms in the first exponent arise from the integration over the quasipar- 
ticle degrees of freedom in eq. (p2D. They represent the density of energy and entropy 
of the quasiparticle gas. 



^ ^qp ^ = ^E' + 36^/,^] , (26) 

St = ^E' [4' In /I -/I In /I + Kft^^ft-ft^^ft)] ■ (27) 



The thermal occupation probabilities are of Bose-Einstein form,[] fl = {e^k/^ _ 
with fl = 1 + fl, and similarly for Furthermore, e| and denote the frequencies 

^This result appears when the field configurations considered are restricted to those for which 
any given quasiparticle mode contains an integral number of elementary quanta, each having the 
energy ek, thereby reducing the functional integral over the fluctuations to a multiple sum over 
occupation numbers ■ This condition is appropriate in the present context where the underlying 
system is of quantal nature; treating the occupation numbers as continuous variables would lead to 
a Maxwell-Boltzmann form, fk = exp(— Cfc/T). 
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of a given quasiparticle mode, = /c^ + /x^, where /i^ is the appropriate eigenvalue of 
the quasiparticle mass tensor M, 

/ij = A(30g + 3 ^ ^ +3^ 5(1)1 y -f^) , (28) 
= A(02 + ^ 50j ^ +5^S4>ly -v^) . (29) 

In these expressions the fluctuations are given self-consistently, 



E'# = E'i^TT^- . (31) 



k 



In the last two expressions in eq. (|25|) , the exponent has been rearranged for 
instructive purposes. The quantity Ko{ipo) = |'?/'q is the bare kinetic energy density 
associated with the evolution of the order parameter. Furthermore, 

VV(0o,Xo) = V^o + + ^ ^qp ^ , (32) 

can be regarded as an effective potential for the order parameter. Finally, the entropy 
density St{4>o) expresses the effective number of accessible quasiparticle states for a 
given magnitude of the order parameter (at the specified temperature T), and the 
free energy density associated with a given order parameter cp is then = Vt — TSt- 
Combining the above results, we note that this latter quantity is given by 

Ft{±) = ^{<Pl-v'f - i/0ocosxo 



^A 50j +2^ y^5(t)ly + 5 ^ 501 y^] (33) 



+ ^ ^'[ln(l - e-'l/^) + 3 ln(l - e"^'^/^)] . 

k 



It is then elementary to calculate the statistical weight ( psD as a function of the order 
parameter, yielding the function Wt{iPo, (l>Oy Xo)- 



3 Dynamical scenarios 

We describe in this section how field configurations corresponding to various instruc- 
tive physical scenarios can be prepared and treated within the linear a model. As we 
have above, we assume throughout that the system is enclosed in a three-dimensional 
torus, i.e. a rectangular box with periodic boundary conditions. 

3.1 Sampling 

We first briefly summarize how it is possible to utilize the above mean-field treat- 
ment of equilibrium to devise a simple approximate method for preparing chiral field 
configurations that represent chiral matter in statistical equilibrium 
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The first task is to sample the order parameter {if), (p) on the basis of the statistical 
weight Wxi'ipo, (po, Xo) given in eq. (^). This quantity factorizes, due to the additive 
form of the exponent. The time derivative ipo is then governed by a four- dimensional 
normal distribution, ~ exp{—QiljQ/2T). Furthermore, since the distribution 

of the magnitude 0o can be pretabulated with the aid of the explicit expression ( p3D 
(ignoring the H term), the associated sampling task is computationally simple. Once 
00 has been picked, the disorientation angle Xo is easy to sample from its conditional 
distribution, Pxixo) ~ exp(— if^o cos Xo)? and the 0(3) direction {■^q,(Pq) of tt is 
uniform on An. 

Once the magnitude of the order parameter is known, the thermal quasiparticle 
distributions are fully determined and the number of quanta in each elementary mode 
is readily sampled, using P{n\^ ~ exp(— nke/t/T) for each of the four principal chiral 
directions. Since the quasi-particle mass tensor is aligned with the 0(4) direction of 
the order parameter, {xo^^o^'^o)^ ^ subsequent 0(4) rotation of (t){r) and li^ir) must 
then be performed in order to express the field in the usual (a, tt) reference system. 

By employing this simple sampling method, it is thus possible to sample field con- 
figurations corresponding to macroscopically uniform matter at a given temperature. 

3.2 Finite sources 

The above method is appropriate for the study of uniform matter. When more 
realistic scenarios are considered, the hot system of physical interest is embedded in 
the normal vacuum in which the order parameter is (f)^^^ = (/^, 0) and -0^^^^ = (0, 0). 
Such systems can be equally easily prepared by applying a simple modulation of the 
sampled fields. 

Let the shape of the hot part of the system, the "source" , be described by the mod- 
ulation function g{r) which drops from unity well inside the source to zero outside. 
For the simple case of a spherical source of radius R we may use 

where a governs the surface width. More complicated source configurations can be 
described similarly. 

To obtain a suitable configuration {il){r),(l){r)) describing the finite source, we 
first pick a field configuration ('0'(r), 4>'{r)) corresponding to uniform matter at the 
desired temperature, for example by using the method described above. This field is 
then subjected to a simple spatial modulation, 

0W = 0vac + ^7(r)[0'(r)-0_] , (35) 
V'(r-) = V'vac + ^7(r)[V'(r)-V'vac] • (36) 

The resulting a field configuration now describes hot matter at the specified temper- 
ature when analyzed inside the source region and it smoothly approaches the vacuum 
value outside. It can then be used as an initial condition for either the full equation 
of motion (^ or its mean- field approximation (^0) . 
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3.3 Pseudo-expansion 

In order for disoriented chiral condensates to appear, it is necessary to allow the 
system to undergo a non-equilibrium relaxation towards the vacuum configuration. 
In high-energy collisions, such a development is expected to occur as a consequence 
of the cooling generated by expansion and radiation. 

Various expansion scenarios have already been considered in the literature. Wang 
and Huang |^ sought to emulate the standard longitudinal Bjorken expansion by 
employing a boost-invariant grid in the beam direction while ignoring the transverse 
dimensions. The natural variables are then the local proper time, r = yi^ — z^, and 
the spatial analog of the rapidity, ?7 = | ln[(t + z)/(t — z)]; the evolution of the field, 
(f){ri, r), was obtained by solving the correspondingly transformed equation of motion 
using random initial conditions for the field strength at each grid point. Later on, the 
treatment of the longitudinal expansion scenario was modified by Asakawa et al. [[1^ 



who included the transverse variables {x,y), while requiring the field to be strictly 
boost invariant, (pix, y,T). A more complete treatment of the longitudinal expansion 
scenario was developed by Cooper et al. who took full account of the transverse 



degrees of freedom in addition to the Bjorken expansion, thus obtaining 0(a;, y, rj, r). 
Analogous expansions in three dimensions have also been considered. Gavin and 
Miiller derived an approximate equation of motion for the order parameter in such 
a scenario by employing a simple time dependence of the decaying field fiuctuations. 
A full inclusion of the expanding three-dimensional geometry was made recently by 
Lampert et al. to leading order in the 1/N expansion, yielding 4>{r], 6, 0, r) with 
V = ^ + r)/{t — r)] being the radial rapidity. 

Since realistic expansion scenarios are both computationally cumbersome and 
physically intransparent, we employ in the present study the simple pseudo-expansion 
method suggested in ref. |jl8|. It consists of augmenting the equation of motion (|l]) by 



a Rayleigh dissipation term of the form —{D/t)dt4> {^^ ~Dil^]f^/t in the k representa- 
tion) and starting the dynamical evolution at to = 1 fni/ c (as is commonly done). This 
term acts as a cooling agency by causing the field fiuctuations to decay in the course of 
time. The associated decrease of the energy density is given hj E = —{D/t) < ip"^ >. 
Moreover, at sufficiently large times, the quasiparticle amplitudes fall off as ~ t~^^'^, 
as can easily be seen by exploiting the approximate decoupling that occurs when the 
field fiuctuations become small. Consequently, the quasiparticle number density (and 
the associated power density) decreases as ~ t~^ as t oo, as is characteristic of 
an expansion in D dimensions. In order to understand that the cooling term can be 
regarded as generating a pseudo-expansion, it is helpful to note that temporal part of 
the d'Alambert operator, d^, transforms into d"^ + {D/T)dr for a Bjorken-type expan- 
sion in D dimensions. Finally we note that the Rayleigh cooling term is unaffected 
by the Hartree factorization and so the pseudo-expansion can be included into the 
mean-field dynamics in the same manner. 
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4 Observables 



After the above developments, we now turn to the calculation of observables relevant 
to the experimental detection of disoriented chiral condensates. We concentrate on 
the pion power spectrum, the pion correlation function, and the distribution of the 
neutral pion fraction which are those that have received the largest attention so far. 

We consider a cubic box with a side length L = 5 — 10 fm and sample the initial 
field configurations from a thermal ensemble at the temperature Tq = 400 MeV, as 



briefly described in sect. |3.1| . The L values are of typical nuclear magnitude and so 
the considered system may be regarded as roughly approximating the interior of the 
excited region formed in a high-energy collision. Moreover, based on the studies in 



ref. |T^, the initial temperature Tq is chosen so as to enhance the chance for the 
system to enter the unstable region; lower values will lead to more stable evolutions. 

After its preparation, the system is subjected to a pseudo-expansion in D dimen- 
sions. The associated attenuation of the field fluctuations then causes the system 
to experience a non-equilibrium relaxation towards the normal vacuum. When the 
fluctuations have subsided sufficiently, and the order parameter has come close to its 
vacuum value, their further exponential evolution can be obtained analytically. It is 
then straightforward to extrapolate to t — >^ cxd and thus extract observable quantities. 

4.1 Power spectrum 

We consider first the power distribution of the pions. Generally, invoking the devel- 
oped mean-field formulation, the total power density at any point in time during the 
dynamical evolution is simple to express, 

Ptotai = ^[V'k ° V'k + 0k ° 0k + 0k ° ° 0k] , (37) 

where M should be replaced by M_ for the term with k = 0. From an observational 
perspective the interest is focussed on the asymptotic stage of the evolution, t oo. 
At that late stage, the field deviates only little from its vacuum value and the different 
quasiparticle modes are effectively decoupled. For the pion field we then have 

7r(r,t) = E^ke^'^-^-*-'^*, (38) 

k 

where the frequency is given by the free pion dispersion relation, = + m^. The 
spectral distribution of the power density for the emerging free pions is then given by 

dPt, 1 

-j^ = ■^Y.^kK-'^kS{uJk-m^- E) , (39) 

where E is the specified kinetic energy of the observed pion. In thermal equilibrium 
the corresponding quantity would be the power density associated with the modes 
perpendicular to the order parameter, 

J p 1 

^ = ^ h S{e, -f^^-E), (40) 
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relying again on the mean-field approximation. 

The results are summarized in fig. ^ The upper-left panel shows the equilibrium 
power density (^) for the initial temperature of Tq = 400 MeV. The solid histogram 
is for the employed side length of L = 5 fm, while the dashed histogram is for an 
eight times larger cube, L = 10 fm. The slight dependence on the magnitude of the 
volume is primarily due to the shell effects associated with the smaller cube. 

The remaining three panels in fig. |T| show the power spectra ( |55D resulting from 
pseudo-expansions in various dimensions, D=l,2,3. For each value of D, 100 initial 
field configurations have been sampled from the thermal ensemble having a temper- 



ature of To = 400 MeV. As already observed in ref. |T^, for relatively slow cooling 
rates, D^l, the evolutions are nearly adiabatic and so there is little effect on the power 
spectrum. When the dimensionality of the expansion is sufficiently large, D^2, the 
system is effectively quenched into an unstable configuration. The low-energy pion 
modes then experience an exponential growth leading to a corresponding enhance- 
ment in the power spectrum |^, [l^, |T3|, [l^ . (This low-energy enhancement causes 
an overall reduction elsewhere, since the displayed spectra have been given a common 
normalization, for convenience.) The resulting distortion of the power spectrum has 
been advocated as a possible signal for DCC formation . 

The results obtained for _D= 1,2,3 with the mean-field dynamics, eqs. (p|-pl]), are 
included as the dashed histograms in fig. |l[ Generally, these results are quite similar to 
those obtained with the unapproximated equation of motion (|^) and so the dynamical 
distortion of the relaxing chiral field appears to be well approximated by the mean- 
field equations. 

4.2 Correlation function 

It is particularly interesting to calculate the correlation function of the pion field. 



since this quantity determines the spectral distribution of the emerging pions |21 



26| , |27| , |28| . We shall therefore extract the equal-time pion correlation function at 
large times when there is no longer any time dependence, except for the overall decay 
caused by the cooling term, 

C^(ri2) = <57r{r + ri2,t)-57r{r,t)> / <6n{r,t)-S7z{r,t)> . (41) 

Here < ■ > indicates a spatial average over the position r within the volume of the 
torus, Q, as well as a subsequent average over the ensemble of histories considered. 
The normalization is such that the correlation function is unity at zero separation, 
C^(0) = 0. The translational symmetry of the scenario implies that the spatial 
dependence is via the separation ri2 = Ti — T2 only. Moreover, to the extent that 
there is invariance under spatial rotations, only the magnitude ri2 = |r"i2| enters. In 
principle, the rotational symmetry is broken by the cartesian lattice, but this effect 
is relatively unimportant. 

For each field configuration considered, the correlation function has been extracted 
in two different ways. One way is to calculate C^^^ = (C^ + C^ + C^)/3, where C!^{ri2) 
has been obtained by assuming that the relative separation ri2 is aligned with the x 
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axis, and analogously for and C^. The other way yields which is extracted 
by aligning with the (111) diagonal. For a cube with side length L, C^^^ has a 
periodicity of L (thus effectively limiting the separation ri2 to L/2), while C^^^^ has a 
periodicity that is \/3 times larger, thereby making better use of the limited volume 
of the lattice. The difference between the two extracted correlation functions may 
serve as an indication of the combined effect of the incomplete rotational symmetry 
arising from the cartesian lattice and the finite number of histories considered. 

The resulting correlation function for the transverse quasi-particle modes in ther- 
mal equilibrium is shown in the upper-left panel of fig. 0, for a number of box sizes, 
L = 5, 8, 10 fm and using a temperature of Tq = 400 MeV. The correlation function 
is neither sensitive to the box size nor to the particular method of extraction, thus 
giving confidence in the method employed. Moreover, it is in good agreement with 
the analytical result for infinite matter (see ref. WW )- 

The three other panels of fig. ^ show the correlation functions extracted at large 
times from the same systems as those considered for the above calculation of the power 
spectrum. Since, roughly speaking, the correlation function and the power spectrum 
are related by a Fourier transformation, the preceding discussion already suggests 
what behavior to expect. Thus, the correlation function for D = 1 is rather similar 
to the thermal form, in accordance with the small distortion exhibited by the power 
spectrum. As the dimensionality of the pseudo-expansion is increased, the weight of 
the low-energy modes grows and a noticeable tail emerges in the correlation function. 
Since the pion modes having the lowest energy are amplified the most, this effect is 
somewhat suppressed by the relative smallness of the box (the lowest quasi-particle 
momentum is given by k^m = L/2tt). This is clearly borne out by the results for a 
larger box, L = 8 fm, which display a more pronounced long-range enhancement. 

It is interesting to note that the mean-field dynamics (solid curves) leads to a very 
good reproduction of the corresponding results obtained with the unapproximated 
equation of motion (solid symbols). This constitutes additional support for the utility 
of the approximation. 

The correlation functions can be conveniently characterized by their full width 
at half maximum. For L=8 fm, this quantity is extracted as 1.0 (1.0), 1.45 (1.30), 
2.75 (2.60) using the exact (approximate) equations of motion with D = 1,2,3, 
respectively. These values reflect the sensitivity of the DCC phenomenon to the 
cooling rate. 



4.3 Fraction: idealized sources 

It has been proposed 0, ^, ^, Q that the occurrence of disoriented chiral condensates 
may be revealed by observing the neutral pion fraction, / = no/n, where rio is the 
number of neutral pions observed and n = n_ + no + n+ is the total pion number. 
In the extreme idealized case where the pion field is uniquely oriented in a particu- 
lar 0(3) direction {i.e. the pion field has that particular orientation throughout the 
source region), only pion modes with that particular isospin alignment can be pro- 
duced. Since that particular direction of the pion field is entirely random, due to the 
isospin symmetry, it is then a matter of elementary geometry to derive the resulting 
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probability density for / 0, ^, 



m) = ^ . (42) 

This result is shown in the first panel of fig. |^. It deviates qualitatively from the 
binomial form expected for a normal scenario in which the charge states of the emitted 
pions are determined essentially independently. Consequently, the distribution of the 
neutral pion fraction, P{f), may provide a suitable signal for the DCC phenomenon. 
However, when assessing the prospects for a measurable signal to actually appear, it 
is important to consider less idealized scenarios, as we shall do below. 

The simple inverse square-root behavior (^2|) occurs only when a single domain is 
formed within which the pion field is fully aligned. However, as a first step towards 
more realistic scenarios, is instructive to consider several different domains, each one 
having its own particular isospin orientation. We shall first consider such sources, 
assuming for simplicity that each contributes the same total number of pions, n. 

When only two domains are considered, N = 2, the resulting probability density 
is easy to calculate by means of an elementary convolution, P2 = Pi * Pi, 

P,{f) = 2 J dfdf" Pi{f) Piif) 5if' + f"-2f) = I I _ 2^ 1^ ^ fj , (43) 



where cos 6 = 1/ \/2f . It is also shown in fig. ^. By performing successive 
convolutions, P/v = -P/v-i * -Pi, it is then straightforward to obtain P(/) for larger 
values of N . The resulting probability densities are shown in fig. |^ for up to A^=6 
sources of equal size, beyond which the trend is clear and fairly uninteresting. The 
mean value of the neutral fraction remains equal to one third for all values of A^, while 
its variance decreases steadily, cr| = 4/45A^. The anomalous behavior for a small 
number of sources quickly reverts to a Poisson-like distribution as several sources 



contribute. These curves are the same as those recently presented in ref. p9[. We 
also note that the calculations made for fig. 3 in ref. for two and three idealized 
sources agree well with the exact results (taking account of the Monte Carlo error). 

It is interesting to see how the above result is modified when the emission has 
a more statistical character. For that purpose we assume that the pion multiplicity 
distribution has an exponential form, for each of the three charge states and for each 
of the iV sources independently. One might regard this modification as a simple 
means of modeling an ensemble of sources having different sizes (as measured 
by the multiplicity of pions contributed). Since the overall normalization of the 
multiplicity is immaterial for the present analysis, the (relative) multiplicity v can 
then be considered as a continuous variable. Thus, the normalized ttq multiplicity 
distribution from a single source would be Pi(z^o) = exp(— i/q) which has unit mean 
and variance. The normalized multiplicity distribution arising from similar sources 
can again be obtained by successive convolution leading to Poisson distributions of 
ever higher order, p7v+i(t'o) = i^o^exp(— z/o)/A^!. In a given event, the neutral fraction 
is then f = Vq/v where = z/_ + z/q + is the total number of pions emitted from the 
A^ sources in one particular event. It is elementary to show that the corresponding 
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distribution is P[{f) = 2 — 2/. Furthermore, for two sources, 



P,'(/) = (P{*PO(/) = { |(F_/)^3^'^^'^ (/>!) ' (44) 

and increasingly complicated algebraic expressions hold for larger A^. The mean value 
remains equal to one third, of course, and the associated variance is aj = 1/18A^. 

The resulting distributions Pj^if) are included in fig. ^ (dashed curves). There 
is a remarkable similarity between the distributions obtained in the two different 
scenarios, although P'{f) has a width that is smaller than that of P{f) by a factor 
of y^5/8 ~ 0.79, consistent with the more statistical nature of the emission. Other 
idealized emission scenarios have been considered. In particular, Cohen et al. 
considered the effect of quantum-mechanical isospin correlations, and very recently 



Anselm and Ryskin calculated P{f) for a rapidly varying source field. The 
various analyses lead to the general expectation that the neutral fraction can display 
an observable signal only if it arises from very few independent sources. 

4.4 Fraction: dynamical simulations 

After the above analysis of idealized cases, we consider now neutral fractions that have 
been generated by dynamical simulations of the type leading to the power spectra 
and correlation functions discussed above, namely pseudo-expansions of thermally 
sampled chiral fields in one, two, and three dimensions. The results are based on 
samples of 100 configurations in a cubic torus with L = 5 fm, the same as those 
employed for the calculation of the power spectra and the correlation functions. 

The neutral pion fraction is sensitive to event-by-event fluctuations in the indi- 
vidual multiplicities and a proper treatment must therefore invoke the underlying 
quantum state represented by the given field configuration. The simplest manner of 
making such an association is to assume that the given pion field 7r(r,t) represents 
a coherent many-body state |3^, |7r > ~ exp(X]k''''k ■ >■ This form leads to a 



Poisson multiplicity distribution and, as a consequence, it contains arbitrarily large 
charge fluctuations, which appears to be somewhat unrealistic. More complicated 
associations displaying more realistic features have also been suggested (see ref. p7| , 
for example) and they lead to various types of correlated emission. An instructive 
analysis of several suggested associations was given recently by Amado and Kogan 



|2^ . In the present work, we assume for simplicity that the multiplicities are given 
by the corresponding expectation values, so that the number of neutral pions con- 
tributed by a particular mode k is ri^ = < > ~ C(jfc|7r^^''p and the corresponding 
total number of pions arising from that mode is rzk = < "n-k > ~ ^^fcTTk ■ T^k- The 
above analysis of idealized sources, sect. E^, suggests that the inclusion of incoherent 



multiplicity fluctuations would tend to degrade the DCC signal. 

Let us first discuss the results obtained with the unapproximated equation of 
motion (p. The calculated distributions of the neutral fraction, P{f) are displayed 
in fig. ^ (left panels). When extracting / we have considered different classes of pions, 
corresponding to various idealized experimental cuts. The most restrictive cut admits 



14 



only those pions that have zero momentum {k = 0), which arise from the isospin 
oscillations of the order parameter, 7r(t). The resulting probability density for the 
neutral fraction (solid dots connected by solid lines) follows very well (to within the 
numerical sampling errors) the anomalous l/v7 form corresponding to the idealized 
one-source scenario (thin solid curve) . This feature is a simple consequence of isospin 
symmetry and the fact that only the single mode having vanishing k contributes. 

The second class of pions consists of those having a kinetic energy E below 200 
MeV. This inclusion of the lowest-lying modes with finite momenta already signifi- 
cantly modifies the distribution of the neutral fraction and it now exhibits a broad 
peak with a maximum near one third. Since the number of pions arising from the 
oscillations of the order parameter is relatively small, as illustrated by the power 
spectra in fig. |l], they have little influence on the resulting neutral fraction, once even 
the lowest quasi-particle modes are included as well. This feature is illustrated by the 
similarity between the dashed and the dotted curves which were obtained by either 
including or excluding the zero-momentum pions, respectively. Although the result- 
ing distribution is still anomalous due to its large width, it no longer displays the 
striking divergence near / ~ and so it does not present a qualitative signal. 

The third class of pions considered includes all the produced pions without regard 
to their energy. The associated neutral fraction distributions are relatively narrowly 
peaked, and more so the slower the cooling is. They are centered around one third and 
the widths resemble what would be expected for statistical production of the pions. 
For example, the distribution obtained for D=3 resembles what would be obtained 
with the statistical emission model considered above, using about nine similar sources. 
We also note that the present results are in accordance with earlier studies p| in 
which a simplified dynamical equation was used to demonstrate that large systems 
[L ^ 16 fm) produce binomial distributions whereas sufficiently small systems (L=4 
fm) display neutral fractions similar to those obtained for ^ 3 idealized sources. 

When considering the dependence of P{f) on the dimensionality of the pseudo- 
expansion, we observe that the anomalous behavior grows more pronounced as D 
increases. This is a consequence of the fact that the anomalous multiplicity distri- 
butions are primarily generated by the low-energy modes and they experience an 
increasing degree of amplification as the cooling rate is increased |T^. Even the dis- 
tribution based on all the pions is noticeably wider for the fastest expansion scenario. 

In order to assess the utility of the mean-field approximation as a tool for dy- 
namical studies, we have treated the same three pseudo-expansion scenarios with the 
approximate dynamics and the resulting / distributions are included in fig. ^ (right 
panels). The agreement between the two sets of results is remarkable. In fact, there 
appears to be correspondence within the errors arising from the finite sample size (100 
evolving field configurations have been considered in each scenario). Therefore, for 
the extraction of this type of observable, the mean-field approximation may indeed 
be useful. 

It is perhaps most striking that the results for all the different expansion scenar- 
ios are so similar. This feature suggests that the appearance of the neutral fraction 
distribution is largely decoupled from the profile of the power spectrum. Indeed, our 
simulations suggest that the distribution neutral pion fraction, P{f), may exhibit 
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an anomalous behavior even if the expansion is farily slow. However, as soon as 
pions with finite momenta are included, P{f) has the appearance of a broad peak 
rather than a decreasing function, so it is qualitatively similar to the normal behavior 
and careful quantitative analysis is then called for in order to identify the DCC phe- 
nomenon. It should also be kept in mind that since the anomalous behavior of P{f) 
is driven by the low-energy modes, the practical observability of the phenomenon 
depends on the abundance of low-energy pions. This inherent feature poses several 
problems. One is the simple geometrical fact that there are relatively few pion modes 
having small momenta, so the signal-carrying pions are relatively scarce. Another 
problem is that such pions are relatively difficult to detect because they move close 
to the beam. Moreover, especially in collider experiments, they may be below detec- 
tion threshold. A rough indication of the relative multiplicity of low-energy pions is 
provided by the power spectra in fig. |l| and, consequently, the amplification produced 
by fast expansions are still expected to be helpful, perhaps even essential for the de- 
tectability of the DCC phenomenon. The effective dimensionality of the expansion 
(or more directly: the rate at which the field fiuctuations subside) is thus a question 
of key importance to the prospects for observing disoriented chiral condensates in the 
laboratory. 



4.4.1 Equalization 

We present here a further analysis that serves to better illuminate why it is so rela- 
tively hard to produce real DCC "domains" during the dynamical evolution. 

We first note that the initial field configurations tend to possess a large degree 
of isospin symmetry. This is especially evident when they are sampled in thermal 
equilibrium at high temperature, but it is expected to be generally the case, since there 
are yet no isospin-breaking mechanisms in play at that early time (if there were, then 
they should better be included in the DCC studies). The only deviations from perfect 
rotational invariance in isospace is associated with the statistical fiuctuations in the 
population of the individual quasi-particle modes. It then follows that the 0(3) tensor 
describing the fluctuations perpendicular to the order parameter, -< 6(f)j_6(f)^ y, is 
approximately isotropic. Consequently, the same is true for the perpendicular part 
of the quasi-particle mass tensor M. 

It is then important to understand how such initial imperfections develop in the 
course of time, as the system expands and cools. It is only the lowest pion modes 
that may be amplified in the course of the pseudo-expansion [|1^], as is borne out by 



the resulting power spectra displayed in fig. Since the contribution of the lowest 
modes to the field fiuctuations is relatively small (even after any amplification they 
may have experienced), the feedback of such amplification onto the mass tensor is 
generally small. The mass tensor can then be regarded as approximately isotropic 
throughout and the different isospin components are then amplified about equally. 
There is thus no tendency for "domain growth" , i.e. a catastrophic isospin alignment 
of different quasi-pion modes. In fact, any isospin imbalance in the mass tensor will 
quickly be counteracted by the amplification (see below): If a given isospin direction 
has larger-than-average fiuctuations, the associated effective mass squared is then also 
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larger than average, and the amphfication rate is correspondingly reduced, thereby 
bringing the fluctuation tensor closer to isotropy. Thus, any initial anisotropy will 
at best be preserved, but is generally degraded by the dynamical effect of the self- 
interaction. 

In order to further elucidate this equalizing mechanism, one may analyze an ex- 
treme (and unrealistic) scenario in which the dominant pion strength is carried by 
unstable modes. It is then elementary to show that any anisotropy in the fluctuations 
would quickly decay. To appreciate this important feature, it suffices to consider a 
simpler world in which the pion field is an 0(2) vector, tt = (7r_, 7r_|_). In this schematic 
model, let the initial fiuctuations be given by < 6it'^ I < ^''^'^ > (1 =1= e), where 
the small parameter e measures the anisotropy in the isoplane. The corresponding 
effective quasi-particle masses /i± are then given by 

jjl = A[< (7^ > + < dn^ > (2 ± e) - v"^] = fil ± eX < 671^ > , (45) 

and the dispersion relations are u"^ = ul ± eX < 671"^ >, where ul = k"^ + /ip with 
Ho being the effective mass for e=0. Inside the region of instability (where u"^ < 0) 
the fluctuations are then amphfied as < 671 > ~ exp(2|u;-|-|t). Consequently, the 
anisotropy decays approximately exponentially. 



The associated characteristic rate is ^ 2A < Stc"^ > /\uJk\ for a given magnitude 
of the wave number. Inserting numbers that are typical of realistic scenarios, we find 
that the decay time is only a fraction of one fm/c. Since this time is relatively short 
in comparison with the characteristic duration of the amplification process (which is 
1-2 fm/c ll^), the dynamics will tend to suppress rather than enhance any anisotropy 
in the isospin distribution of the quasi-pions. 



5 Concluding discussion 

The possible formation of disoriented chiral condensates in high-energy hadron and 
heavy-ion collisions has generated considerable interest in the past few years, since 
the associated anomalous pion multiplicity distributions may offer a means for testing 
our understanding of chiral symmetry. In order to elucidate the conditions for the 
occurrence of DCC phenomena and the prospects for their experimental detection, 
it is necessary to carry out dynamical simulations of the non-equilibrium evolution 
experienced by the chiral field as it relaxes from an initially very excited state, in 
which chiral symmetry is approximately restored, towards the normal vaccuum in 
which the symmetry is significantly broken. 

The most popular tool for such dynamical studies has been the linear a model in 
which the chiral degrees of freedom are described by an 0(4) real classical field subject 
to a simple non-linear interaction. Even though this description is relatively simple, 
ignoring all other degrees of freedom (such as those represented by other meson type 
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or individual quarks), it still presents a significant computational challenge. It is 
therefore of practical interest to develop useful approximate methods for solving the 
non-linear equation of motion for the chiral fields. 

In the present paper we have presented a simple mean-field approximation that 
is based on the Hartree factorization technique and retains all the terms in a 
expansion. It describes the system in terms of an evolving order parameter and a gas 
of independent quasi-particles whose interaction is via their self-consistent medium- 
modified effective mass tensor. This can be regarded as a dynamical extension of 
the treatment developed in ref. |^ in the context of statistical equilibrium and it 



presents a simple and useful conceptual framework for analyzing the complicated 
DCC dynamics, as already illustrated in ref. |jl8||. Because of the effective decoupling. 



the mean-field dynamics can be solved very efficiently in the k representation, leading 
to a savings in computational effort by over an order of magnitude.^ 

The present treatment thus relies on a succession of approximations which can be 
briefly summarized as follows. The general framework is the familiar linear a model 
in its classical form where real a and tt fields experience an interaction containing up 
to the fourth power of the fields. Furthermore, the fields are subjected to periodic 
boundary conditions and their initial form is chosen to reflect Bose-Einstein statistics 
(see below for a further discussion). On top of these general features, the present 
work develops a "mean-field" approximation in which the quasi-particle degrees of 
freedom interact only via their self-consistent effective masses. 

It is instructive to note the analogy of the present treatment with the semi-classical 
mean-field models that have proven very successful in the context of nuclear reaction 
dynamics at lower energies. The closest analogue is the Vlasov approximation in 
which the individual nucleons move independently in a common self-consistent mean 
field (as the motion of the quasi-particles in the present treatment is governed by 
the self-consistent mass tensor). In order to take account of their underlying quantal 
nature, the initial phase space distribution is specified in accordance with the Fermi- 
Dirac statistics governing real nucleons (as the present treatment employs the Bose- 
Einstein statistics appropriate to the quasi-particle meson gas). This type of approach 
is justified by the fact that the time scales associated with the inevitable reversion 
to Maxwell-Boltzmann statistics are long in comparison with the typical durations of 
the processes under study (whether nuclear reactions or DCC formation). 

In addition to presenting the formal development of the mean- field treatment, we 
have also sought to illustrate its practical utility by testing its accuracy for a num- 
ber of observables of particular interest in the quest for observing DCC phenomena. 
The comparisons with the corresponding results of the unapproximated non-linear 



^For the unapproximated equation of motion the non-linear character of the interaction renders 
the k representation inefficient and the x representation is preferable, even though the Fourier 
transformations associated with the initialization and the spectral analyses are relatively costly. The 
computational effort is then roughly proportional to the number of lattice points, (L/c?x)^, where the 
lattice spacing is dx ~ 0.2 fm. In the mean-field approximation, there is no direct coupling between 
the different modes and the k representation is then the economical choice. The associated effort is 
proportional to the number of modes included, {2K + 1)"^, where L/K « 1 fm is the minimum wave 
length considered. The resulting gain factor is then ~ (5/2)"^ and exceeds ten. 
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equation of motion suggest that the developed mean-field approximation is indeed 
quantitatively useful for this kind of simulation. 

The calculated results are interesting in their own right, as they illustrate various 
key features of the DCC phenomenon. We employed a cubic volume with a side 
length of L = 5 fm as representative of source sizes expected in heavy-ion collisions. 
A larger source would behave similarly in the present idealized kind of calculation. 
In particular, since the DCC anomaly is associated primarily with the low-energy 
pions {i.e. the signal- carrying pions have close to zero momentum in the frame of 
the source), the ideal 1/a/7 behavior of P{f) will appear only if the observation 
is restricted to A; = 0. Consequently, as the box is enlarged this single component 
will constitute an ever smaller fraction of the total power. This feature illustrates 
how the collection of pions from a large spatial volume tends to wash out the DCC 
signal. A cleaner signal emerges for a smaller volume, so a key question is how small 
phase-space domains can be experimentally separated. 

Moreover, we have concentrated on evolutions prepared at Tq = 400 MeV, in order 
to increase the chance for the system to enter the unstable regime where exponential 
amplification occurs |T^. If the initial temperature were lower, the magnification of 
the low-lying pion modes would be smaller as well, or entirely absent, but the results 
for the pion fraction would not change appreciably, as is borne out by the insensitivity 
of P{f) to the value of D. 

Once the system has cooled sufficiently, the different modes effectively decouple. 
Each mode k contributes pions having an independent orientation in isospace, TTk. 
The resulting pion multiplicity distribution can then be understood as resulting from 
a collection of such independent sources, in a scenario similar to that discussed in 



sect. [4.3| . Accordingly, the resulting probability density P{f) is the corresponding 
convolution of all the contributions from all the individual modes, with proper ac- 
count taken of their different "sizes" as measured by the contribution to the total 
multiplicity. As an instructive aid to the interpretation of the dynamically calculated 
distributions of the neutral pion fraction, P{f), we extended various earlier discus- 
sions of idealized emission scenarios, considering a number of similar but independent 
pion sources each contributing either a fixed or a stochastic multiplicity. This analysis 
brings out the difficulty in obtaining the ideal l/vT signal in realistic scenarios when 
there are several different sources of the observed pions. 

Simulations including more realistic features are clearly called for. One important 
limitation of simulations done with the standard linear a model is the absence of the 
many other degrees of freedom in the system. It is of interest to understand their 
effect on the DCC phenomenon. Moreover, the excited systems produced in high- 
energy collisions are embedded in the normal vacuum and they are born in a state 
of rapid expansion. Since this latter feature may be essential for the occurrence of 



DCC phenomena ||T^, it appears desirable to incorporate realistic initial expansion 
patterns. These must be obtained on the basis of more fundamental descriptions, such 
as parton dynamics, and until such information has become available the simulations 
of DCC phenomena can only be suggestive. 
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Figure 1: Pion power spectrum. 
The upper-left panel shows the power spectrum of the transverse quasi-particle modes 
for thermal equilibrium at the temperature Tq = 400 MeV. The system is enclosed in 
a cubic torus with size L = 5 fm (solid histogram); for comparison is also shown the 
results for L = 10 fm (dashed) (the irregularities are caused by the shell structure of 
the cube). The other three panels show the power spectrum of the pions emerging 
asymptotically as a result of propagating 100 field configurations that were sampled 
from the thermal distribution at Tq. The solid histograms were obtained by employing 
the exact equation of motion (|l|) while the dashed histograms are for the mean-field 
approximation, eqs. (|9|-|TTD, in addition to the Rayleigh cooling term —{D/t)dt<p which 
roughly emulates a Bjorken-like expansion of matter in D dimensions (ordinary static 
equilibrium is thus recovered for D=0). The normalization is such that the total 
power is the same in all cases. The power carried by mode having k = (the order 
parameter) is indicated by the solid circle and square for the exact and approximate 
dynamics, respectively. 
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Figure 2: Pion correlation function. 
The pion correlation function C7r('^i2) defined in eq. PT|). The upper-left panel is the 
correlation function for the transverse quasi-particle modes in thermal equilibrium at 
the temperature Tq = 400 MeV, obtained with various values of the box size, L=5,8,10 
fm. The other three panels were obtained at large times after 100 field configurations 
were sampled from that thermal ensemble and subjected to a pseudo-expansion in 
D=l,2,3 dimensions, using either the exact equation of motion (|l]) (solid curves) or 
the mean- field approximation, eqs. (|9|-|lTl) (dashed curves). The results obtained with 
the exact dynamics for a box with L=5 fm are shown by the solid symbols, as obtained 
by directing the relative separation ri2 either along the lattice (squares) or diagonally 
(circles). Exact results obtained with L=8 fm are indicated by the corresponding open 
symbols. The results of the corresponding approximate dynamics are shown by the 
solid and dashed curves for L=5 fm and L=8 fm, respectively. 
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Figure 3: Neutral fraction from idealized sources. 
The normalized probability density, P{f), for the neutral fraction / = Uq/u obtained 
for N similar sources, assuming either the exact same multiplicity contribution n from 
each source (solid curves) or an exponential distribution of the individual multiplicities 
n_, no, n+ for each source separately (dashed curves). 
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Figure 4: Neutral fraction from dynamical simulations. 
The normalized probability density, P{f), for the neutral fraction f = rio/n obtained 
by subjecting a thermal sample of 100 field configurations to a pseudo-expansion in 
D dimensions, using either the exact equation of motion (|l]) (left) or the mean-field 
approximation (right). The various curves were obtained by considering different 
classes of pions when extracting /: those pions having k = only (dots connected by 
solid lines), those with kinetic energies E below 200 MeV either excluding or including 
k = (dots connected by dotted or dashed lines, respectively), and all of the pions 
(dotted curve). For reference is shown the result for the idealized case of a single 
source (light solid curve), also shown in the upper-left panel in fig. ^. 
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